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Abstract 

A new version for the relativistic generalization of the wide class of static 
quark-ant iquark confining potentials is suggested. The comparison of this ap- 
proach with other ones, known in literature, is considered. With the use of 
Logunov-Tavkhelidze quasipotential approach the first-order retardation correc- 
tions to the heavy quarkonia mass spectrum are calculated. As expected, these 
corrections turn out to be small for all low-lying heavy meson states. 
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1 Introduction 



Within the framework of the constituent quark model it is quite natural to study 
the properties of bound qq systems on the basis of Bethe-Salpeter equation. De- 
spite the remarkable success in the quantitative description of the bound-state 
masses and formfactors, achieved with the use of instantaneous (static) inter- 
action kernel in this equation [1-4], the completely relativistic approach to the 
problem is still lacking. From the physical point of view one can expect that the 
dynamical retardation corrections (i.e. the corrections, coming from the explicit 
dependence of the quark-quark interaction kernel on the relative energy variables) 
to the bound-state characteristics must be small for heavy quarkonia and may 
become significant in the light-quark sector. Moreover, one expects smooth static 
limit, when the masses of the constituent quarks tend to the infinity. In prac- 
tice, however, the situation is more complicated, owing to the infrared-singular 
behaviour of the (phenomenological) confining quark-quark interaction kernels. 
Since at the present stage the exact derivation of such relativistic kernels di- 
rectly from QCD is not known, the different prescriptions are assumed for the 
ad hoc relativistic generalization [5-11] of the static interquark confining poten- 
tials. As to the one-gluon exchange part of the interquark potential, which has 
proven to be significant in the quantitative description of meson data, it can be 
uniquely generalized to 4 dimensions with the use of the field-theoretical argu- 
ments. In difference with the one-gluon exchange potential, the "relativization" 
of the confining potentials, in general, introduces a new mass parameter in the 
theory [0, [|, [10], which may be fixed, using the additional constraints either 
directly on the relativistic counterpart of this potential 0, || [K], [Il| or on the 
bound-state equation ||. Neither of these constraints can be preferred from the 
physical point of view, rendering the identification of the dynamical retardation 
effect ambiguous. 

The existence of the additional free parameter in the theory stems from the 
necessity of the infrared regularization of the " confining" kernel in 4 dimensions. 
As a result, the smooth static limit, in general, does not exist || |Tl[ and the 



dynamical retardation corrections to the bound-state characteristics turn out to 
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be large in the heavy-quark sector || [TT|, rendering doubtful even the concept 
of the confining interaction in 4 dimensions ||. Moreover, the straightforward 
relativistic generalizations for some widely used confining potentials (e.g. for 
the oscillator potential) seem do not possess the discrete spectrum |12|] in the 
entire energy region |lTJ. Thus, the further investigation of the problem in the 
framework of Bethe-Salpeter equation with the relativistic kernels is needed in 
order to clarify the notion of the confining force in 4 dimensions. In other words, 
one has to answer the question, whether there exist relativistic kernels which are 
compatible with confinement, lead to a smooth static limit and, therefore, can be 
expected to be obtained from QCD with the use of nonperturbative methods. 

2 The Relativistic Generalization of the Static 
Confining Potentials 

In the present work we suggest a new prescription for the relativistic general- 
ization of the wide class of static confining potentials, where some difficulties 
which are inherent for the prescriptions, suggested earlier [5-10], are presum- 
ably avoided. In order to clarify the physical meaning of this prescription, let us 
consider the case of two space-time dimensions. We start from the power-law po- 
tentials V a (r) = r a . The corresponding Bethe-Salpeter kernel in two dimensions 
has the form: 

K a (x) = K a (x , Xl ) = 8(x )V a (\xx\) = 8(nx)(-x 2 ) a/2 (1) 

where and n = (1, 0). We assume that the relativistic generalization 

of this expression is the " Lorentz-average" over all possible directions of the unit 
vector n which can be obtained from the vector (1, 0) by the proper Lorentz 
transformations. In 2 space-time dimensions this vector transforms as n' ((p) = 
hq cosh ip — ri\ sinh ip, n'^ip) = n\ cosh ip — no sinh ip, ip = - and — oo < <p < oo 
and the group-invariant measure is dip. Then the "average" over all possible 
directions of n can be defined as follows: 

K a (x) = {-x 2 r' 2 < S(nx) > n = (-x 2 ) a/2 jj /_ oo d<p6{x ri (<p) - x^ip)) (2) 
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where M is an arbitrary Lorentz-invariant normalization constant. The integral 
in the r.h.s. of eq. (2) can be easily evaluated: 

/oo 
dip5(x n' (ip) - xm'M) = 9(-x 2 )(-x 2 )- 1/2 (3) 
-oo 

The constant M is fixed from the condition that in the static limit the kernel 
(2) reduces to its static counterpart V a (r): j'^dx^K^XQ^Xi) = V a (\xi\). This 
gives: 

K a (x) = l ) a J (-^)(-^) la - 1)/2 (4) 

y/nL (2 + 2) 

Thus, the relativistic kernel (4) in two dimensions can be interpreted merely 
as the static potential, averaged over all possible directions of the quantization 
axe, which are equally acceptable from the physical point of view. For the case of 
4 space-time dimensions the straightforward evaluation of the integral over the 
group-invariant measure in analogy with the eqs. (2)-(3) is not possible. Despite 
this the "average" of the static kernel (1) over all possible directions of the unit 
vector n can be defined for this case as well, and we again arrive at the expression 
(4) (for the details see Appendix A). Note that the Coulombic interaction in 4 
dimensions (a — — 1) is excluded in the expression (4). 

Having obtained the expression (4), we can forget about its "derivation" and 
consider (4) merely as an anzats for the relativistic generalization of the power- 
law potentials V a (r) = r a . Further, we consider the exponential static potential 
Vexp{ r ) — e _Atr . Expanding this potential in powers of r and using eq. (4) for 
each term of this sum, we obtain the prescription: 

e-» r - e(-x 2 ) (-^v^) + i(-* 2 )- 1/2 1^2(1; \ \; -^-)) (5) 

where J u and p F q denote, respectively, the Bessel and hypergeometric functions. 
Eq. (5) enables one to apply the procedure of the relativistic generalization to 
the wide class of potentials which can be written in the following form: 

V(r) = J dfiC{fi)e-^ r (6) 

where C(/x) must obey to the certain conditions in order to provide the conver- 
gence of the integral over d/i in the relativistic case. 
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Now we discuss the particular case of linear confinement (a = 1) in detail. 
Note that the anzats similar to (4) for the gluon propagator was employed earlier 



in ref. [fT"3| , |23|. The Fourier transform of eq. (4) for a = 1 reads: 



Kl{q) ^ Wo " (151 + iO)') 3 + (gg - (\q\ - iO)*y) (7) 
i.e., in other words, in order to obtain the kernel (7) in the ladder approximation, 
the principal-value prescription must be used in the gluon propagator instead of 
the familiar causal one. One can find here a close analogy with the 2-dimensional 
QCD (see, e.g. [14-16]), where the gluon exchange is definitely known to confine 
quarks as well as with the gluon propagator from ref. |2T| . As it is well known |23 



+ ig 2 C F j,S(p)T^p,p) J -^D(q) 



due to the principal-value prescription in (7) the most severe infrared singularities 
in the equations for the Green's functions are avoided. To demonstrate this, we 
consider the Schwinger-Dyson equation for the quark propagator written in the 
Feynman gauge: 

S-\p) = Sv\p) + ig 2 C F J t^jT^Cp - V) T "(P ~ ^P)D(q) = 

= S \p) + ig 2 C F J {S{p - q)T»(p - q,p) - S(p)T»(p,p)) D(q) + 

d A q 

where S(p) and So(p) are full and free fermion propagators, r M (fc,p) is the dressed 
quark-gluon vertex function, D(q) is the dressed gluon propagator, g is the quark- 
gluon coupling constant and Cp is the quadratic Casimir operator in the funda- 
mental representation. If the principal-value anzats (7) is used for D(q), then the 
last term in eq. (8) is finite and we obtain the subtracted form of the Schwinger- 
Dyson equation, where the infrared singularity is softened. In difference with this, 
in the Euclidean formulation of Schwinger-Dyson equations, which implies the 
causal prescription for all propagators after the continuation to the Minkowski 
space, the last term in eq. (8) diverges and requires the infrared regularization 
. Consequently, the relativistic generalization of the static interquark interac- 
tion kernels introduces an additional mass parameter when the causal prescription 
is used in these kernels or, equivalently, when one works in the Euclidean space 



from the beginning. However, no such parameter is needed, when one uses the 
principal-value prescription (7) or, in general, the prescriptions (4)-(6). In the 
absence of the additional free mass parameter one expects the existence of the 
smooth static limit for the calculated meson observables, while the corrections to 
the static limit should allow for the unambiguous evaluation. 



3 First— Order Quasipotential for the qq Systems 

The formulated approach for the construction of the relativistic confining kernels 
below we apply to the calculation of the dynamical retardation corrections to 
the bound qq system masses. The static confining potential was taken to be 
linear + constant term: V c (r) = kr + c and the spin structure was chosen to be 
the equal-weight mixture of scalar and the fourth component of vector: O c = 
® I2 + 7° <8> 7°) which is perhaps the simplest choice from the more general 
ones [1-4] and provides the existence of the stable discrete energy levels. The one- 
gluon exchange part of the potential was neglected since we were interested in the 
retardation corrections, coming from the confining part of the potential. Further, 
for the simplicity, we have restricted ourselves to the equal-mass case. We have 
used the Logunov-Tavkhelidze quasipotential approach |H| in order to reduce 



the original, 4-dimensional Bethe-Salpeter equation to the 3-dimensional one, 
which can be numerically solved with the use of the conventional mathematical 
methods. 

The quasipotential equation for the 3-dimensional equal-time wave function 
(pip) (in c.m.f.) is written in the following form [|T0| , [11]: 

[M B - hip) - h 2 i-p)Mp) = -* 7l °7 2 ° \ J pHM fl ;p,#(g) (9) 

where Mb is the mass of the bound qq system and hi = ctiPi + rwyf, m being the 
mass of the constituent quark. Keeping in mind that the retardation corrections 
for heavy quarkonia are expected to be small, further we restrict ourselves to the 
first-order quasipotential formalism, where V_ is given by the following expression 



1C, 11 



V (1) iM B ;p,q)=<p\G 'GoKGoQv 1 \q> (10) 
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Here Go is the free two-fermion Green's function, K is the Bethe-Salpeter 
equation kernel and the procedure A for any operator A is defined as follows: 

A(P;M= [%*A(P;p, q )p. (11) 

and 

Go = G o7 ?7 2 n; U = (A«A« - A^A^h?^; 

A& = B±!±. Wi = (m*+$)i (12) 
The partial-wave decomposition of the eq. (9) is done with the use of the 



standard technique ]2|, [Hj, PL1| . In the wave function (pip) we retain only the 
contribution from the "double-positive" component tfi( ++ \p) = A^ A^ <f(p) , as 
the role of the negative-energy component is shown to be small provided the 
stable solutions exist 0. Using the substitution 



2w J \ ai P I V 2w J \ - a2 P 
and decomposing the Pauli spinor x (p) partial waves 

X (+ \p) = E <ft\LSJMj>R^j(p); n = V -- 5 = 0,1 (14) 

LSJMj P 

from the eq. (9) we obtain 

4 f°° o , , 4 A m 2 \ 



[Mb - 2w(p)}R[t(p) = -J™ q 2 dq [V^ip, q)-U 



w(p)w(q) 



(15) 

where, for the case of local quasipotential 

V[ +) (P, q) = J r 2 drj L (pr)V c (r; M B )j L (qr) (16) 

]l being the spherical Bessel function and we, as in ref. |[10|, TT[] , have neglected 
small mixing between L = J — 1 and L = J + 1 states. 

After doing some algebra with the projection operators the positive-energy 
component of the first-order quasipotential, corresponding to the interaction ker- 
nel (4) with the use of eqs. (10), (11) can be written in the following form: 

r ( i + — 1 

V {+) (M B ;p,q) = f d z xe~ l{ P~^ ^ — ^r Q x 



Pi 

x / dr(l — r 2 )^ a ~ 1 '^ 2 {9{T)e i ^ MB ~ w ^~ w ^ rT + 9(— T )e~ i ( MB ~ w(j ''~ w ^' )rT ) (17) 

Neglecting relativistic corrections in the exponentials Mb — w(p) — w(q) = 
Mb — 2m + 0(— ) = — e#, as well as the imaginary part of this expression in 
analogy with the refs. [[IT], |TT|, |T§||, we arrive at the local first-order quasipotential: 

V {+) (r;e B ) = T V ^ N r a C dr(l - t 2 )^' 2 cos e B rr = 

■ v^rQ + f) J-i 



2r 


a/2 






r 











Hence, the relativistic generalization of the static potential V" c (r) = kr + c 
gives the first-order local quasipotential 

fc(r, e fl ) = + cJoCear) (19) 

which accounts for the retardation effect and reduces to V c (r) in the limit — * 0. 

As it can be seen from eq. (19), the account for the dynamical retardation 
effect in the case of pure linear potential effectively leads to the colour screen- 
ing at an intermediate distances. It should be pointed out that such behaviour 
qualitatively agrees with the results of calculations for the unquenched lattice 
fermions in QCD [p0[| . At a larger distances the deviation of the retarded poten- 
tial from the static one becomes significant, and one can no further rely on the 
first -order calculations. Note, however, that for the case of heavy quarkonia the 
latter difficulty causes no trouble since the wave function of the qq bound system 
in this case rapidly vanishes with the increase of r and, therefore, does not "feel" 
the oscillating "tail" of the potential at a large distances. 



4 Results and Discussion 

It is not obvious from the beginning whether the potential (19) leads to the 
discrete energy levels due to its oscillating behaviour at r — > oo. Let us, there- 
fore consider the equation (15) with the potential (19) in detail. Passing to the 
nonrelativistic limit and neglecting for a moment the "constant" term in (19), 
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proportional to c, in the configuration space we obtain the following differential 
equation: 



f"(z) + (a cos 2z + &)/(*) = (20) 

where f(r) = rR(r), R(r) = Ro(r) being the radial wave function of the bound 
state in the configuration space (for simplicity we assume the angular momentum, 
L = 0),z= i((M B -2m)r-7r/2), a = k(M »J m) , b = and the boundary 

conditions imposed on f(z) are / (— f) = and /(+oo) = 0. 

The equation (20) has been extensively studied in the mathematical physics 
(see, e.g. [p2fl ). We shall remind some results of this investigation. Namely, if 
fi(z) is the particular solution of the eq. (20) with the following initial conditions: 

A(0) = 1, f[(0) = (21) 

and 

cosh27r/i = /i(7r) 
Then the general solution of the equation (20) has the form: 



/(*) 



de^ip^z) + C 2 e- 2 » z ip 2 (z); cosh27r/i > 1 
{Ci cos 2^2 + C 2 sm2vz)ipi(z) + (C 2 cos2i/^ — C\ sm2uz)ip2( z )', 

| cosh27r/i| < 1; fi = \v 
^ Ctf?P*tp 1 (z) + C 2 e- 2pz tp 2 (z); cosh 2tt/i < -1; n = p + \ 



(22) 

ip\(z) and (p 2 (z) being the periodic functions in z with the period n. 

Due to the fact that we consider the equation (20) on the semi-infinite inter- 
val — f < z < +oo, it is possible to find the normalizable solutions decreasing 
exponentially at z — >• +oo (Ci = and | cosh27Tyu| > 1, eq. (22)). The eigenvalue 
condition then reads as: 



; r2 



0, | cosh27Tyu| > 1 



(23) 



Thus, the equation (20), despite the oscillating behaviour of the potential at 
the spatial infinity, allows for the discrete spectrum provided |cosh27r/i| > 1, 
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corresponding to the condition Mb — 2m < in the limit \Mg — 2m\ << 2m. 
Adding the constant term, proportional to c, it is natural to suppose that, for a 
small \Mb — 2m\ the discrete energy levels exist for Mb — 2m — |c < 0. Thus the 
potential (19) in the nonrelativistic limit acts like the potential well. Note that 
the similar potential (the rising potential screened at large distances, r > 1 Fm 
was successfully used for the description of meson spectrum in the framework 



of coupled Dyson-Schwinger and Bethe-Salpeter equations e.g. in ref. p4 
Therefore we expect that the equation (15) with the potential (19) gives the 
reasonable description of the low-lying meson states. 

At the next step we have attempted to solve the eq. (15) numerically, expand- 
ing the unknown radial wave function R^sjip) i n the complete orthonormalized 
basis of the nonrelativistic oscillator wave functions 0, |H], [LJJ 



R^AP) = Po 3/2 E cttsjRndp/Po) (24) 



>< + ' ,„\ - ,,- ; V2^(+i 
n=0 

where 



_L 



(25) 

and po is an arbitrary scale parameter. Inserting (24) in the equation (15) and 
truncating the sum at some fixed value N max , we arrive at a system of linear 
algebraic equations for the coefficients c^ SJ . If the procedure converges with 
the increase of N max , the eigenvalues Mb are determined from this system of 
equations. The calculations show that the final results do not depend on the 
scale parameter po, but the appropriate choice of this parameter leads to the 
faster convergence of the series (24). It should be stressed that if the solution of 
the equation (15) does not exist, this reveals in the divergence of the procedure 
with the increase of N max despite the fact that the potential matrix elements are 
calculated in the exponentially damping wave function basis. 

Since the potential (19) depends on the unknown binding energy, €b = 2m — 
M B , of the qq system, the equation (15) is solved with the use of the iteration 
method. Namely, we solve the equation with the static potential, V c (r) = kr + c 
and determine the eigenvalues Mg^- At the next step these static values are 
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substituted into the potential (19) in order to determine the corrected spectrum 
which, in its turn, is used, as an input, in the next iteration. We have checked 
that, typically after 10-15 steps, the iteration procedure converges for the most 
low-lying heavy quarkonia energy levels. 

In the table 1 the results of calculations of the dynamical retardation correc- 
tions to the heavy quarkonia mass spectrum are presented. In these calculations 
the parameters k and c were taken to be k = 0.21 GeV 2 , c = —1.0 GeV. The 
constituent quark masses were chosen to be m c = 1.72 GeV and rrib = 5.10 GeV 
in order to fit J/ip and T masses. As we see from table 1, this set of parame- 
ters gives the reasonable description of heavy meson mass spectrum in the static 
approximation. As expected, the dynamical retardation corrections turn out to 
be small (typically a few percent) for all low-lying quarkonia states given in this 
table. 

As we see, the present approach to the relativistic generalization of the static 
confining potential does not suffer from some of difficulties inherent of the conven- 
tional ones [5-10]. Namely, an additional mass scale parameter does not appear 
in the "relativized" counterpart of the static potential. Moreover, this approach 
possesses a smooth static limit and yields small corrections due to the retardation 
in the heavy quarkonia mass spectrum. This, unlike the results of refs. [f|, |TT| , 
does not manifestly contradict with our expectation at least in the heavy-quark 
sector. An important remark, however, should be made, concerning the approach 
suggested in the present paper. Namely, the kernels off the type considered here 
(eqs. (4), (5)) can never be obtained in the Euclidean formulation of QCD (see eq. 
(7)). Consequently, such fundamental properties, as microcausality and unitarity, 
do not necessarily continue to hold in this approach. Of course, one can appeal 
again to the 2-dimensional QCD, where, in spite of the principal-value prescrip- 
tion used in the gluon propagator, the uncoloured current-current correlation 
functions are known to have the correct analytical structure | 16|| . Note, however, 
that in 2 dimensions the gluon field is not a dynamical one [15], unlike the case of 
4 dimensions, where the different components of the gluon field (dynamical and 
nondynamical) may be subject to the different boundary conditions. Here one 
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finds an analogy with the result obtained in the ref. |2T| where, in particular, 
it was demonstrated that the ghost degrees of freedom (in the covariant gauge) 
can be eliminated from the theory at the cost of Feunman's boundary conditions 
in the propagator for timelike and longitudinal gluons, being replaced by the 
principal-value prescription. An alternative point of view was presented in ref. 
[23|1 , where the quantization procedure for the gauge field is carried out with the 
account of the fact that this field is never observed in in— or out— states. As a 
result, the principal-value prescription emerges in the propagator for all compo- 
nents of the gauge field instead of the conventional Feunman one. Consequently, 
the anzats (4), corresponding to the "minimal" relativistic generalization of the 
static interquark interactions, can be regarded as an appropriate candidate for 
the truly confining interquark interaction kernel in 4 dimensions. 

The authors are thankful to Dr. B.Magradze for many helpful discussions. 
One of the authors (A.R.) acknowledges the financial support under grant UR- 
94-6.7-2042. 

Appendix A 

Below we define the "average" of the static kernel (1) over all possible direc- 
tions of the unit vector n in the 4-dimensional Minkowski space. This approach 
can be applied to the case of any space-time dimensions and, in particular, to 
the case of 2 dimensions, considered in the text. 

We start from the Lorentz-invariant regularization of the "average" of 5- 
function, entering the expression (1): 

— <e^ nx) - eX \9(x 2 ) + 9(-x 2 ))> n =5 + + 5- (A.l) 

-oo lH 

Let us consider first the case x 2 > 

5 + = 6(x 2 )(x 2 )- 1/2 \im f°° 7^< e a(m?+) ~ eA2 > n , rft = x^/Vx^; rf + = 1 (A.2) 

Expanding e lX( - nv+ ^ in powers of A and defining the "average" in a familiar 
manner: 

<%!-•• >n= for k odd 
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< 1 > n = 1; < n m n M2 > n = -g^ 2 ; etc (A.3) 

we obtain 

< e lA ^+) > n = 1 + + • • • = F+(A) (A.4) 

In order to calculate F + (X) we consider the 4-dimensional Euclidean space 
and define the average: 

< e **<»V> >„ £ = 1 / dV^(n s2 - l)e iA ^); n^ 2 = ijf = 1 (A.5) 

Using the relations 

< nf ■ ■ - nf. >„e= for fc odd 

< 1 > n E= 1; < njnj 2 >„ B = ^<W 2 ; etc (A.6) 

it is easy to ensure that (A.5) coincides with F + (X). Passing in the integral in 
(A.5) to the angular variables and substituting the result into (A. 2), we finally 
obtain 

5 + = 0(x 2 )(x 2 )- 1/2 ^ o ^ J°° d\ J 1 dyy/l - y 2 e iXy - tX2 = 9(x 2 )(x 2 )-^ 2 ^ 

(A.7) 

The case x 2 < can be considered in a similar way. However, due to the 
change of sign rf_ = —1, where rft = x^/ V—x 2 , in the integrand e tXy is replaced 
by e Xy 

5_ = e{-x 2 ){-x 2 )- l ^ 2 \im^— I" dX f 1 dyJl - y 2 e Xy - eX2 (A.8) 

e^O AlX 2 J-oo J -I V 

and, as a result, £_ diverges in the limit e — > 0. Including this divergent (Lorentz- 
invariant) factor in the overall normalization constant Af, we see, that only the 
term <5_ survives in the limit e — * 0. Finally, substituting (A.8) into the eq. (2) 
and determining the normalization constant in the static limit, we again arrive 
at the expression (4) for the relativistic Bethe-Salpeter kernel. 
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Mesons 




N 2S+1 Lj 


1 


2 


3 


4 


cc 

r? c (2.980) 

^(3.590) 

J/^(3.097) 

^'(3.686) 

/i c i(3.526) 

Xco(3.414) 

Xci(3.511) 

X c2 (3.556) 

bb 

Vb 

V'b 

T(9.460) 
T'(10.023) 

T'"(10.355) 

hi 

X6o(9.860) 
Xfti (9-892) 
X6 2 (9.913) 

Xw(10.232) 
Xw(10-255) 
X6 2 (10.268) 


o-+ 
o-+ 
\— 
1— 
1+- 

0++ 

1++ 

2++ 

o-+ 
o-+ 
1— 
1— 
1— 
1— 
1— 
1+- 

0++ 

1++ 

2++ 
1+- 
0++ 

1++ 

2++ 


I 1 S 
2% 

2 3 S' 1 
l^i 
l 3 Po 
l 3 Pi 

1 3 P 2 

1% 

2% 

l 3 Si 

2 3 S 1 
l 3 D l 
3 3 S 1 
2 3 D 1 
l l Pi 
l 3 Po 
l 3 Pi 
1 3 P 2 
2 1 P 1 

2 3 P 
2 3 Pi 
2 3 P 2 


3.095 
3.690 
3.096 
3.691 
3.445 
3.445 
3.445 
3.445 

9.463 

9.899 

9.463 

9.899 

9.938 

10.250 

10.277 

9.720 

9.720 

9.720 

9.720 

10.097 

10.097 

10.097 

10.097 


3.204 
3.739 
3.204 
3.742 
3.460 
3.460 
3.460 
3.461 

9.619 

9.966 

9.619 

9.966 

9.993 

10.255 

10.291 

9.835 

9.835 

9.835 

9.835 

10.109 

10.109 

10.109 

10.109 


0.109 
0.049 
0.108 
0.051 
0.014 
0.014 
0.014 
0.014 

0.156 
0.067 
0.156 
0.067 
0.055 
0.005 
0.014 
0.115 
0.115 
0.115 
0.115 
0.013 
0.013 
0.013 
0.013 


3.5-10~ 2 
1.3-10- 2 
3.5-10- 2 

1.3- 10 2 
4.2-10 -3 
4.2-10" 3 
4.2-10 -3 
4.2-10~ 3 

1.7- 10- 2 

6.8- 10 -3 

1.7- 10- 2 

6.8- 10- 3 
5.5-10- 3 

5.1- 10 -4 

1.4- 10 -3 

1.2- 10- 2 
1.2-10 -2 
1.2-10 -2 

1.2- 10 -2 

1.3- 10 -3 
1.3-10 -3 
1.3-10 -3 
1.3-10- 3 



Table 1. The dynamical retardation corrections to the heavy quarkonia mass 
spectrum. 1) The meson mass in the static approximation, Mg*' (GeV), 2) The 
meson mass with an account of the retardation effect, M% e€) (Gev), 3) The size 

of the retardation correction, M% et) - (GeV), 4) | M% e£) - M^l/M^ . 
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